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Abstract 

A two-component model is developed for diffractive electroproduction of p°, <f> and Jflf, 
based on non-perturbative and perturbative two-gluon exchange. This provides a com- 
mon kinematical structure for non-perturbative and perturbative effects, and allows the 
role of the vector-meson vertex functions to be explored independently of the production 
dynamics. A good global description of the vector-meson data is obtained. 
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1 Introduction 



High-energy exclusive photo- and electroproduction of vector mesons offers a variety of insights 
into the diffractive mechanism. The choice of different vector mesons and a range of pho- 
ton virtualities allows one to move from the primarily nonperturbative regime to the primarily 
perturbative within one framework, and to explore kinematical regions where neither is predom- 
inant. Vector-meson production also has the benefit of a high rate, but it has the disadvantage 
of dependence on the choice of the vertex functions which couple the vector mesons to the qq 
pairs. In both respects it differs from deep virtual Compton scattering, which is theoretically 
better defined but has the disadvantage of a small cross section. 

On the basis of the factorisation theorem 0, exclusive vector meson production can be con- 
sidered as three separate processes: the fluctuation of the (virtual) photon into a qq pair; the 
interaction of the qq pair with the proton; and the formation of the vector meson from the 
qq pair which naturally involves the vector-meson vertex function. It may be argued that the 
structure of the nonpertubative vector-meson vertex function invalidates the proof of the fac- 
torisation theorem as it leads to additional contributions. However it has been shown ||, at 
least in a simple model for the vertex function, that gauge invariance ensures that the additional 
contributions cancel and factorisation is preserved. 

The aim of this paper is twofold: to obtain a global description of exclusive vector meson photo- 
and electroproduction; and to explore the choice of vertex function. 

The interaction of the qq pair with the pomeron is modelled by two-gluon exchange, figure 1, 
which can be applied both to nonperturbative and perturbative gluon exchange. The former is 
based on the model of Diehl J| and the latter either by utilising the gluon structure function 
or following the model of Cudell and Royon Q . 

The approach has the advantage of providing a common kinematical structure in which it is 
possible to separate the part of the vector-meson production amplitude describing the kinemat- 
ics from the part which describes the dynamics of the process ||. This separation then allows 
the nonpertubative and perturbative contributions to be combined at the amplitude level. Thus 
the approach follows recent ideas about two pomerons and two-component models of diffraction 
through |12j| , combining "soft" (nonperturbative) and "hard" (perturbative) terms. This is 



essential for a global approach as neither a nonperturbative nor a perturbative model alone can 
describe all the observed features of diffractive vector-meson photo- and electroproduction. 

In a previous paper |5j we have successfully applied this approach to calculate p electroproduc- 
tion in a model which avoids the vertex function complication, following |T3|| , by considering 
open gg-pair production with the invariant mass of the qq pair restricted to the region of the p 
mass. However this method cannot be applied to higher-mass vector mesons: for a qq with an 
invariant mass in the region of the p there are no states available except the p itself (making 
a 10% allowance for ui production), but this is not true for the and J/ijj. To obtain a global 
description of p, <fi and J/ip photo- and electroproduction a kinematical framework involving 
vector-meson vertex functions has to be used. 
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Figure 1: Two of the four vector meson production diagrams in the kinematical framework of 
Cudell and Royen [|J. The other two diagrams differ by reversal of the quark charge flow and 
give the same contribution to the cross section. The off-shell quarks are marked by crosses, 
dashed lines indicate cuts along which the quark lines are put on-shell, the minus sign indicates 
the momentum of an antiparticle. 



In figure 1 the quark lines marked with a cross are off-shell, and the vertex function must 
necessarily take this into account. We use the prescription of Cudell and Royon M for these 
vertex functions. They appear very different from the usual on-shell vertex functions, for 



example those of Brodsky and Lepage ||14|| . However they are both obtained by a boost of the 
same vertex function from the vector-meson centre-of-mass and differ only by the imposition 
or relaxation of the on-shell condition. 

Kinematical and dynamical considerations of the two-gluon exchange model are coverd in sec- 
tions 2 and 3 and the vector-meson vertex functions are discussed in section 4. The model is 
applied to the data in section 5, and conclusions are given in section 6. 



2 The Kinematical Framework 

The kinematical framework of our gluon exchange models is depicted in figure [TJ and was 
developed by Cudell and Royen As noted earlier, it involves the vector-meson vertex 
function $ (£) with either the quark or the anti-quark off-shell. The various four momenta used 
in our discussion are also defined in figure [l]. In terms of them, the photon virtuality Q 2 = —q 2 
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and the squared centre of mass energy of the photon-proton pair W 2 = (p + q) 2 . The traces 
corresponding to the gg-loop in diagrams (a) and (b) of figure [1] are given by 



XT = Tr{$ (£) ( 7 .e) (m, + 7 . + 7" (m, + 7. [q + k + £ - ±V\) 

x Y (m, +J.[q + £-kV]) (7-e) (m q +r[t~ |^]) } 

T b ^ = Tr{$(£)( 7 .e) (m, + 7. + 7" (m, + 7. [g + A; + £ - ±V 

x ( 7 .e) (m, + 7- [fc + * - \V\) 7 M (m 5 + 7- - \v]) } . 



(1) 



(2) 



One of the two quark lines emerging from the j*qq vertex is off-shell with different virtualities 
in each of the diagrams of figure [l], giving 



q + l — W 



k + £-\V 



m„ 



2£.q - £.V - q.V - Q 2 - ml + \M 2 



-ml = k.k + 2k.£- k.V + £.£-£. V - m 2 + \Ml 



(3) 
(4) 



for the denominators of the corresponding propagators. Here we assume the diffractive ampli- 
tudes are completely dominated by their imaginary parts, which are evaluated using the cuts 
shown by the dashed lines in figure 0. The sum of both diagrams for the gg-loop gives 
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where the 5-functions are due to the on-shell conditions along the cuts of the quark and the 
anti-quark lines respectively. The denominator 



2£ 2 - 2m 2 + \M 2 



v 



(6) 



of the propagator for the off-shell quark forming the vector meson is the same for both diagrams 
of figure ^[ In the Regge limit the proton line gives a contribution Ap a p 13 in amplitude and 
the intermediate proton state is cut and its mass neglected, yielding 8 (\p — k] 2 ^j. The rest 
of the diagram, including the gluon propagators and the description of the interaction of the 
pomeron with the proton, is contained within the dynamical part V, which is model-dependent. 
Formally it is the (gauge-dependent) gluon propagator that contracts the indices at PaPp and 
T^y. Practically, the leading contribution in the Regge region comes from g afl g^ v in the gluon 
propagators. Details of this and of the decompositions of the four-vectors in terms of p and 
q for the Regge region, where W is significantly greater than any other scale present, can be 
found in 0, 15]. Here we would only like to note the following points. 

It is convenient to use the light-cone variables Pt = {P + , Pt, P ) in decompositions and in the 
further derivation, where P ± = P° ± P 3 and the two-vector P t lies in the transverse plane, 
defined as the plane perpendicular to the 7*p axis. The variable z is also often used. It is 
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the fraction of the "+" momentum of the photon carried by the quark, so that P+ = zq + , 
implying = (1 — z) q + for the anti-quark. The decomposition of the gluon four- momenta k 
and [k — A] show j|] that the gluon four- momenta are predominantly transverse, \k 2 \ « k 2 . 

For clarity we rewrite equation (21) of Q using our notation with the dynamical part V and 
the light-cone variable z: 
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(7) 



where we have introduced £ t , k t for twice the transverse parts £ t , k t of the four-vectors £, k. For 
(f) and Jfi> electroproduction, f g = — |, | corresponding to the charge of the quark forming the 
vector meson while the linear combination of the uu and dd quark anti-quark pairs forming 
the p -meson gives f g = The expressions for N L and N Tr are given in Q and 



D 
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The first line in D originates from the denominator of the off-shell quark propagator P , equation 
(H); similarly, the second and third lines are proportional to Pb and P a respectively, equations 
(H) and (0). The differences from the expression for D given in M are due to the fact that we 
attribute the gluon propagators to the dynamical part of the pomeron exchange. 



Assuming s-channel helicity conservation the differential cross section is given by 
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(9) 



where the polarisation of the photon beam e expt is a known characteristic of the experiment. 
For HERA, e eX pt ~ 1- For fixed-target experiments, it typically lies in the range 0.5 to 0.9 
depending on the energy and photon virtuality. 

Finally we note that the kinematical expressions in this derivation and in [|J are valid for the 
pair of diagrams depicted in figure [3]. The other two diagrams differ by reversal of the quark 
charge flow and thus have different traces, cut conditions and propagators than ([l|)-@ leading 
to different decompositions of the four-vectors £ and k. However, the only net difference is a 
change of sign in front of £ in all expressions. Because the integral over £ t is two dimensional, 
(0) gives the same answer regardless of the sign in front of £ and the additional diagrams result 
in a factor of 2 in front of the final expression for the amplitude. 
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Figure 2: The two possibile ways in which the two gluons can couple to the valence quarks of 
the proton. In the nonperturbative LN approach as applied by Diehl || it is argued []TB[ that 



the diagrams in which the gluons couple to different quarks are suppressed and subsequently 
can be neglected: see text. 

3 The Dynamics of Pomeron Exchange 

We investigate two "summation models:" these include both "hard" (perturbative) and "soft" (nonperturbative 
components, since neither alone can account for all the data. Before describing these summa- 
tion models, we first summarise the models for the hard and soft terms on which they are 
based. 



3.1 The "soft term" 

For this contribution we use the nonperturbative approach of Diehl || based on earlier work 
by Landshoff and Nachtmann (LN) . In it the gluons are assumed not to interact with each 
other and a nonperturbative gluon propagator 



T-^np ( k ^ 



M, 



Tip 



k 2 



(n - 1) (4 



(10) 



is used with n 



4. The normalisation M np is determined from the condition 



dk 2 



afv np {k 2 



47T 



(11) 



The phenomenological parameters /3q, which describes the effective coupling of the pomeron to 
the proton, and /i are determined from the total pp and pp cross section data and from deep 
inelastic scattering: (3q ~ 2.0 GeV -1 and /xq ~ 11 GeV [17]. For the nonperturbative couplings 



of the gluons to the quarks forming the p a value atg ~ 1 is taken. Of course the precise value 
of aig cannot be strictly specified, and as we shall discuss in section 5 there is some flexibility 
through the interplay with the choice of vertex function for the vector meson. 

Landshoff and Nachtmann [0 have argued that diagrams in which the nonperturbative gluons 
couple to different valence quarks in the proton, as shown in figure 0, are suppressed and can 
be disregarded. Hence only the diagrams where both gluons couple to the same valence quark 



6 



are calculated. Each of the three valence quarks is incorporated into the proton using the Dirac 
form factor Fi p (t), where t = A 2 . The energy dependence of the soft pomeron comes via a 
factor x F aF ^ in the amplitude, where 

x = M v + Q 2 ~ t (12) 



and a P (t) = 1.08 + 0.25 1 is the soft pomeron trajectory The coupling a s at both vertices 
at the proton end is taken at a nonperturbative scale, i.e. ag is used. For the vertices where 
the gluons couple to an off-shell quark line the coupling is taken at a perturbative scale A 2 = 
(fi + m 2 ) (Q 2 + My) /M v which, gued in [Q, is a typical scale for the whole upper part 
of the diagram. Thus for the dynamical part of the nonperturbative approach one has 

V np = F lp (t) x 1 -"*® af a s (A 2 ) V np (-k 2 ) V np (-[k-A] 2 ) . (13) 

This term alone gives an energy dependence that is too flat at the higher values of Q 2 due to 
the soft pomeron intercept. 

3.2 The "hard term" 

Two alternative models, the standard perturbative QCD approach and the Cudell-Royon ap- 
proach, will be considered. 

3.2.1 The "standard" perturbative model 

The perturbative QCD approach can only be calculated at t — 0, and is based on the ideas of 
Martin et al. |T3|] who applied it to p -meson electroproduction. In it the pomeron is modelled 
as a pair of perturbative gluons, using the perturbative gluon propagator V p {k 2 ) = 1/k 2 . The 
gluons are considered as part of the proton so that there are no ^/aZ, couplings for the two 
bottom vertices. In principle the gluon flux can be obtained from the unintegrated gluon 
density f{x Fl k 2 ), which gives the probability of finding a t-channel gluon with the momentum 
squared h\ in the proton. However, a special treatment of the infrared region is required because 
the unintegrated gluon density f(x P ,k 2 ) is undefined as k 2 — > and numerically unavailable 
below some value of k 2 = Qq, which varies with the parton distribution chosen and is typically 
in the region from 0.2 to a few GeV 2 . The linear approximation as suggested in [DJ is used 
to account for the contribution to the integral from the k\ < Qq region. This procedure has 
no direct physical significance. It serves only to provide a continuous integrand and acts as a 
means of normalisation of the perturbative contribution. A simple cutoff at an appropriate Qq 
would be equally effective but somewhat less elegant. 

As no direct physical significance can be attached to the contribution from this infrared part 
of the perturbative term there is not an element of double counting. The separation between 



7 



"perturbative" and "nonperturbative" is given uniquely by the energy dependence of the two 
contributions. An implication of this approach is that the perturbative (hard) term can con- 
tribute at Q 2 = 0, which is a feature of two-component models. 
Thus for the dynamical part of the perturbative approach one has 



- / 7 x 7 -/(W* 2 [*-A] 2 ) 



(14) 



fc 2 [A; - A] 

where f(x P , k 2 ) is related to the gluon distribution g(x P , Q 2 ) by 

2 /"<2 2 dfe 2 2 
x P g(x P ,Q ) = y —^f(x P ,k t ) (15) 

with the inverse 

f(x P ,k 2 ) = k 2 ^ — 2 L, (16) 

This applies at at t = 0, and the experimental slope or some other ansatz must be used to 
compare with the integrated cross section. Here we merely note that this term alone gives an 
energy dependence which is clearly too steep for much of the data. 

3.2.2 The Cudell-Royen (CR) model 

Cudell and Royen propose a somewhat different approach to the "hard" contribution. Again 
perturbative gluon propagators of the form 1/k 2 are used. However, since the cross section is not 
infrared divergent, it is suggested that the divergencies from diagrams where both gluons couple 
to the same quark are cancelled in the infrared limit by the divergencies from the diagrams 
where each gluon couples to a different quark, figure |2|. In terms of form factors, the former 
case is described by the Dirac form factor of the proton i*i(t); while for the latter the form 
factor 

S 2 (k, k -A) = Ft (k 2 + [k - A] 2 + c k. [k - A]) , (17) 

depending on the momenta of both gluons, is applied with c ~ — 1 as suggested by Cudell and 
Nguyen ||19|| . Unlike the "standard" perturbative approach, this approach does not contain any 



energy dependence but describes the i-dependence at a fixed energy. The energy dependence 
has to be introduced by hand, just as it was for the nonperturbative term, via a factor x P ajP °. 
We assume a flat hard-pomeron trajectory a Fo = 1.44, independent of t. Further the overall 
normalization is not uniquely specified as it is not obvious what value of a s should be used for 
coupling perturbative gluons to bound quarks. Cudell and Royon 0, ^(J introduced an effective 
factor 1Z in the cross section with TZal = 0.6 ||, |2(J, a procedure which we adopt here. In this 
way we finally obtain 

V CR = ^} „.(*»)„.([*- A]') ^VS-A]'"^ ^ ' <18) 
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3.3 The Summation Models 



Here we suggest two "summation models" which combine both hard and soft terms at the 
amplitude level in order to obtain a global description of vector meson production data. 

3.3.1 Summation model SI 

This model is based on our earlier work ||, which modelled p electroproduction by "open pair" 
production in the region of the p mass. In particular, a successful description of the data 
was obtained by combining the nonperturbative amplitude of section 3.1 with the perturbative 
amplitude of section 3.2.1, using an empirical slope parameter to describe the t-dependence of 
the latter. In doing so, we exploited several gluon distributions to calculate the perturbative 
contribution using the PDFLIB program libraries |21| for numerical calculations. However, 



it was found that the best fit to the p -meson electroproduction data, and especially to the 
energy dependence of the production cross section, was obtained using the CTEQ4LQ [^2 



gluon distribution. We continue to use the CTEQ4LQ gluon distribution in the present paper. 

In this paper we explicitly incorporate vertex function effects using (0) in order to treat the 
4> and Jfi> as well as the p electroproduction in a common framework. To do this, we again 
need to extend the perturbative amplitude to t ^ 0. Since the proton in the vector meson 
production process remains intact, we suggest describing the t-dependence at the proton end 
by the proton form factor Fi p (t). In this way we arrive at summation model SI: 

V S i{s,t,Q 2 ) = V np {s,t,Q 2 ) + F lp {t)V p {s,Q 2 ) , (19) 

where V np , V p are given by (|13"D, ( |i~4"D respectively. 

3.3.2 Summation model S2 

Another possibility is to model the "hard" component with the CR term - Pc , _r,(^B|). The t- 
dependence is thus automatically provided, but the CR approach does not specify the energy 
dependence. In order to obtain a global description of vector meson production, Regge energy 
dependence corresponding to the "hard pomeron" |7], 23, || is introduced into the CR term 



by hand, in much same way that it was introduced in the nonperturbative term using the soft 
pomeron trajectory. In this way one obtains summation model S2: 

Vs 2 (s,t,Q 2 ) = V np (s,t,Q 2 ) + V CR (s,t,Q 2 ) (20) 

where V np , Vcr are given by (|l~3|),(|l"8|) respectively. 
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4 Vector Meson Wave Functions 



The choice of the vertex function is crucial in vector-meson production models as it determines 
the virtualities dominating the integral over the gg-loop; the overall normalisation and Q 2 de- 
pendence of the cross section; and the longitudinal to transverse ratios. Unfortunately the 
detailed forms of the vertex functions are unknown and only their general analytical properties 
are established from various constraints f24j . Therefore in practice, the chosen vertex functions 
provide an essentially phenomenological description of the valence quark content of the vector 
meson. Here we shall consider possible forms for the vertex functions by starting from phe- 
nomenological wave functions for vector mesons in their rest frame; and then boosting to the 
light-cone taking into account the off-shell nature of the quark line. 



4.1 Wave Functions in Centre-of-Mass Frame 



25, 26, 27, 28, 29, BQ] of vector-meson wave functions is suggested 



The most popular choice 
by long distance physics. This tells us that a hadron at rest can be described to a good 
approximation as a system of constituent quarks moving in a harmonic oscillator potential 
with a Gaussian wave function 



A/bexp 



2pi 



(21) 



where L 2 is the squared 3- momentum of either the quark or anti-quark, pp is the Fermi mo- 
mentum and M is the normalisation. We investigated five alternatives, the details of which can 
be found in E|. The first is the power-law wave function ^5, ^8 




(22) 



with n = 2 in our case. This was found not to be an acceptable choice as it was not possible to 
obtain even a qualitative description of the data, particularly for the longitudinal/transverse 
ratio for the p and <ft which are rather sensitive to the wave function details. 

The four other wave functions are obtained by solving the non-relativistic Schroedinger equation 
with four different potentials ||31|| . The first three of these are: 

• a power-law potential |32[ 



— a\ 



do I — 



0.1 



(23) 



a logarithmic potential 
b\ + 6 2 log(^) 



(24) 
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a Coulomb-plus-linear potential (the Cornell potential) 



+ -2 + C 3 

c 2 



(25) 



where the a^b^Ci are various mo del- dependent parameters. The fourth is the QCD-inspired 
potential of Buchmiiller and Tye [|55]| , which has a rather complicated position-space form. It 
is linear at large distances and quasi-Coulombic at short distances. The deviations from pure 
Coulombic behaviour reproduce the running of the strong coupling constant, and the global 
shape of the potential is essentially determined by two parameters — the QCD scale A and 
the QCD string tension motivated by the light meson data. Non-relativistic wave functions 
are reliable only if the meson and both constituent quarks are heavy compared to the average 
internal momentum, and it is still not clear whether one should use them for the J/\l/. In partic- 
ular Frankfurt, Koepf and Strikmanf^] have used wave functions from various non-relativistic 
potential models to show that the integration region where the quark's transverse momentum 
is larger than the charm mass can contribute up to one third of the gg-loop integral in J ft 
production. For the p and such wave functions can be still be considered as an alternative 
choice, albeit without any firm foundation. 

Finally, all these vector meson wave functions have to be normalised to reproduce the leptonic 



decay width of the meson T\ 



in the vector meson rest frame. Furthermore, the normal- 



isation has to be calculated for one quark leg off-shell to reflect the kinematics of the vector 
meson production models of figure [TJ. The derivation of the normalisation, and its reduction to 
the on-shell case in the appropriate limit, has been given by Cudell and Royen|Q. Their final 
result is 



64 (a em e q ) 



where 



S(|L|) 



in a ill + p 



2L 2 



Q (|L| 



d\L\ 



(26) 



, Ox 

— ^ $ L 2 



(27) 



\Ml-m\ . 



(28) 



and m q is the constituent quark mass. For u, d quarks we take m q = 0.3 GeV, for strange 
quarks m s = 0.45 GeV and for charmed quarks m c = 1.5 GeV. 

The next step is to transform the vector meson centre-of-mass wavef unctions into the vertex 
functions &(z : \£t\) used in (0). This is done by first rewriting the centre-of-mass wavefunctions 
in terms of invariants, which are then re-expressed in terms of the light-cone variables z and 
£ t to obtain the vertex function. We do this first using the Cudell-Royon prescription || for 
the quark or antiquark off-shell, and then show that if both are put on-shell it goes over to the 



Brodsky-Lepage prescription [14[ . 
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4.2 The Cudell-Royon prescription 



In figure |l] the quark and the meson have four- momenta P q — \V + 1 and V respectively. 
In the meson's centre-of-mass frame, P£ = (E q , L) and V* = (My, 0) so that the invariant 



quantity 



becomes M v E q . This gives 



v 



(29) 



and 



L 2 = E\ - P£P 



i.V 
Ml 



1 2 



(30) 



for the squared three momentum of the quark in the centre of mass frame, where no assumption 
about the on-shell or off-shell nature of the quark has been made. Thus one obtains the relation 



CM 



LC 



e.v 

My 



(31) 



between the vector meson wave function $cm expressed in the centre-of-mass variable L 2 and 
$lc expressed in terms of appropriate invariants. This is the Cudell-Royen(CR) prescription 
HI. One finally has to rewrite the four- vector i in terms of the variables z and £ t , using equations 
given by Cudell and RoyenQQ which take account of the fact that the quark is off-shell and 
the anti-quark is on-shell 



\V + 



(32) 



in the two diagrams of figure [l]. The asymmetry between the quark and antiquark results in 
an asymmetry under the transformation z <-» 1 — z, which interchanges the "+" momentum 
carried by the quark and antiquark. However, when calculating vector-meson production, one 
must also take into account diagrams corresponding to those in figure |T] but with the quark 
charge flow reversed. The cut conditions then put the quark (instead of the anti-quark) on-shell, 
giving the reverse of equation (|32|) . Thus the asymmetry present within each pair of diagrams 
is no longer present once the vector-meson production amplitudes from all four diagrams are 
summed. 

As an illustration, the vertex functions $(z, \£t\) corresponding to the Gaussian and logarithmic 
potential vector meson centre-of-mass wave functions (plD are shown in figures |3] and [|, where 
the asymmetry under z «-> 1 — z is clearly seen. One also sees that the vertex functions are 
centred towards the point z = 0, \£ t \ = 0, implying that the most likely configurations are those 
where the off-shell quark carries only a small fraction zq of the meson's longitudinal momentum, 



t In B our l t is denoted by l t 
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while the on-shell anti-quark carries most of it, namely (1 — z)q. However the vertex functions 
never actually reach the point z — 0, \£ t \ = 0, because the region around it is unphysical, 
corresponding to L 2 < 0. The border of this region correspond to L 2 = 0, where the quark and 
anti-quark have no relative momentum in the meson rest frame. In the other case, with the 
anti-quark off-shell and the quark on-shell, the CR prescription gives similar vertex functions, 
but centred towards z = 1 and \£t\ = instead of z = and \£t\ =0 so that the on-shell particle 
again carries most of the meson's longitudinal momentum. 

In the numerical calculations we focussed mainlyQon Gaussian wavefunctions with values of pp 
ranging from 0.1 GeV to 0.6 GeV for p and mesons and from 0.2 GeV to 1.2 GeV for the Jffy. 
The motivation behind these ranges of values is the transverse size of the corresponding vector 
meson; in particular, that for the light vector mesons is similar to that used in the literature 
for the pion [T2^, ^6], E^]. As we shall see, our best results for electroproduction were obtained 
with pf = 0.5,0.6,1.2 GeV for the p, <fi and Jf$> mesons, for which the appropriate values of the 
normalisation constant are Mg — 3.597, 3.226 and 2.905. 



4.3 Relation to the Brodsky-Lepage prescription 



The prescription of Brodsky and Lepage ]TJj, often found in the literature [f26], |36], |57| , connects 



the wave functions in the centre-of-mass frame and the light-cone frame by equating the off- 



shell propagator e = M 2 — (E™=i ki) 



for a particle of mass M whose constituents i 



in the two frames. In the general case the propagator 



, n have masses m; is given by 



M 2 -(E?=i^) E2=iLi = [CM] 

M 2 - E?=i (it + m 2 ) /zi E?=i iu = ELi zt = 1 [LC] 



(33) 



where Lf = (L°,Lj) and Zj, £u are the centre of mass four momenta and light-cone variables 
respectively. For a system of two-particles of equal masses this gives Li = — L 2 . The relation 
L\ = L\ is then used |26| implying that both constituent particles are on-shell. 

According to the Brodsky-Lepage(BL) prescription the propagators in both frames are equated 
yielding 



it + m\ 
Az(l-z] 



(34) 



and 



$cm (L 2 ) 



LC 



it + ™\ 



Az(l 



— m„ 



(35) 



$A fuller discussion of the various possible wavefunctions is given in u\ 
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Alternatively, one could equate the expressions for invariant mass of the gg-pair in the vector 
meson centre-of-mass and light-cone frames 

Ml = 4 (f +m ;)=|±^ (36) 

yielding the same prescription (@). The functions (|35D are often called light-cone wave func- 
tions; and for quarks and antiquarks of equal mass, the prescription is seen to be explicitly 
symmetric under z <->• 1 — z. 

The BL prescription is not applicable in the kinematical framework of figure [I], m which the 
quark (or antiquark) is off-shell. On the other hand the CR prescription (|30|), which allows for 
this, reduces to the BL prescription if both both particles are put on-shell 

~lV-e] 2 = m 2 q [\V + l] 2 = m 2 q (37) 

yielding 

tV = LI = m 2 - \M V (38) 

but not in general 1.1 = or = 0. Substituting ([58]) into (^) gives 

L 2 = \M V - m\ . (39) 

so that the square of the three-momentum of the quark in meson rest frame is fixed by the 
mass of the meson My and by the choice of the mass of the quark m q . However, L 2 is not zero 
provided the m q is not chosen to be exactly half of the meson mass. Identifying the mass of the 
vector meson My with the invariant mass M x of the gg-pair and substituting (|5E| ) into (|39|). 
one again obtains the BL prescription (|34]). 

If the quark mass is exactly half of the meson mass, Eq.(|3"8"D gives 

tV = £1 = (40) 

corresponding to a gg-pair with zero relative momentum in the meson rest frame. The quark 
and the anti-quark share equally the meson's four-momentum. In the light-cone variables it 
reads 

^ (2) (4) ^ (z - I) (41) 
which has also been used as a very basic light-cone wave function in the literature |(], |39|| . 



Figures ^| and | show vertex functions obtained using the Brodsky-Lepage prescription(i.e. 
light-cone wavefunctions), to be compared with the vertex functions of Figures |3] and ^ which 
were obtained using the Cudell-Royon prescription from same centre of mass wavefunctions. As 
can be seen, Brodsky-Lepage wave functions are symmetric under z <-> 1 — z with a maximum 
at z = |, \£t\ = 0. In other words the most probable configuration is where the quark and 
antiquark share equally the vector meson's longitudinal momentum. Further, only one single 
value( "height") of the wave function, determined by the value of L 2 , enters in any given vector 
meson production calculation since the L 2 is fixed by ([39]) once the masses of the quarks and 
the vector meson are fixed. In calculating (^), this can be imposed via a separate condition 
stating the on-shellness of the quarks. Such a condition, reducing the dimensions of integration, 



is indeed present in the vector- meson production models exploiting the BL prescription [20, 38 
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CR.F 



Figure 3: p -meson relativistic Gaussian wave 
function (CR precription) with p F = 0.6 GeV, 
with M = 2.803. 




Figure 5: p -meson relativistic Gaussian wave 
function pp = 0.6 GeV, BL precription. 




CR.F 



Figure 4: p -meson non-relativistic logarith- 
mic potential wave function(CR precription) 
and M = 2.763. 




BL.F 



Figure 6: p -meson non-relativistic logarith- 
mic potential wave function, BL prescription. 
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5 Results 



We have already said that the power-law wave function ( ^2j) is not appropriate, and we do not 
discuss it further. The results obtained using the wave functions obtained from the power-law 
potential (p3|) , the logarithmic potential (|2~4[) , the Coulomb-plus-linear potential (|5|) and the 



Buchmiiller-Tye potential |35[ are almost identical. They predict successfully the longitudi- 
nal/transverse ratio for each of the p (f> and J/9 and the correct shape for da/dQ 2 in each case. 
However it is not possible to obtain simultaneously the correct normalisation of all the cross 
sections. If the normalisation is adjusted to fit the p cross section, say, then the predicted J/9 
cross section is too high. Conversely, if the normalisation is adjusted to the J/9 cross section 
then the predicted p cross section is too low. The problem is that there is no flexibility in the 
wave functions: they are all fixed by the parameters of the potentials. This is not the case for 
the Gaussian wave function (|2l|) for which the parameter pp can be adjusted independently for 
each case. 

It turns out that the Gaussian wave function corresponding to the best choice of pp is very 
close to the wave functions obtained from the power-law, logarithmic, Coulomb-plus-linear and 
the Buchmiiller-Tye potentials for the p and <ft but differs significantly in the case of the J/9. 
This is shown in figures |7|, [5] and from which it is obvious that for the J/9 the model requires 
a much narrower transverse distribution in configuration space than is provided by the wave 
functions obtained from solving the Schroedinger equation for specific potentials. 

For clarity of presentation we show results only for the Gaussian wave function, and comment 
on the principal differences in results obtained with the alternative vertex functions. The pa- 
rameters in our calculation are the pp for each of the p, <fi and J/9, which control the transverse 
size of the vector-meson wave functions; the momentum cut-off when the gluon structure 
function is used in the perturbative term, or the overall normalisation of the perturbative term 
in the Cudell-Royon model. The normalisation of the nonperturbative term is effectively fixed 
by the p photoprodution cross section which it dominates ||. Two sets of results are shown: 
SI refers to non-perturbative plus gluon structure function; S2 refers to non-perturbative plus 
Cudell-Royon. In both cases the same parameters for the two-gluon exchange contributions are 
used for the p, <fi and J/9, so that the only difference allowed for different vector mesons is in 
the value of the relevant pf. 

We start with the p. The results for da/dQ 2 for and ap/ap with pp = 0.5 GeV are shown 
in figures 10, 11 and 12 for SI and in figures |14| and |H| for S2. The result for da/dQ 2 
at W = 75 and for ap/ap, which has little energy dependence, are satisfactory, but that for 
da/dQ 2 at W — 15 clearly falls below the data at all Q 2 . This is not surprising as at this lower 
energy there is a contribution from reggeon exchange which we have not taken into account. 
The results for da/dQ 2 , other than the normalisation, are not strongly dependent on the choice 
of pf. The shape is reproduced satisfactorily for 0.3 < pp < 0.6 GeV. However ap/ap depends 
very strongly on pp, the ratio rising rapidly with decreasing pp, and restricting the choice of 
Pf to the upper end of the range. The S2 result is slightly better overall. The results obtained 
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Figure 7: Comparison of the Gaussian 
wave function (|2lD for the p with pp = 0.5 
(solid line) with the wavefunctions obtained 
from non-relativistic potentials. The dashed 
lines, from top to bottom, correspond to the 
Buchmiiller-Tye, logarithmic, power-law and 
Coulomb plus linear potentials respectively. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



|L| (GeV) 

Figure 8: Comparison of the Gaussian 
wave function ( pi]) for the with pp = 0.6 
(solid line) with the wave functions obtained 
from non-relativistic potentials. The dashed 
lines, from top to bottom, correspond to the 
Buchmiiller-Tye, logarithmic, power-law and 
Coulomb plus linear potentials respectively. 
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|L| (GeV) 

Figure 9: Comparison of the Gaussian wave 
function ([21]) for the Jf$> with p F = 1.2 
(solid line) with the wave functions obtained 
from non-relativistic potentials. The dashed 
lines, from top to bottom, correspond to the 
Buchmiiller-Tye, logarithmic, power-law and 
Coulomb plus linear potentials respectively. 
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using the wave functions derived from the specific potentials are very close to those shown, and 
provide an equally satisfactory description. 

It is clear from figures [H] and [13| that an increase in normalisation would convert a good 
description of the data into an excellent one. If we were considering p meson production in 
isolation then it would be appropriate to do this by increasing otf\ However, as we shall 
see, this would then impact adversely on the cross sections for the and Jfty, particularly 
on the former by making it too large. The problem of simultaneously obtaining the correct 
normalisation for the photo- and electroproduction of the p and the within the constraints of 



pre-defined wave functions is well-known; for example see [40 



The results for the with pp = 0.6 GeV are shown in figures [Tjj, [17| and [18] for SI and in figures 
|1^ , |20| and ^ for S2. The results in both cases are very good. It is particularly satisfying that 
the description of the low-energy data is good, as in this case there is no reggeon contribution. 
Once again the results obtained using the wave functions derived from the specific potentials 
are very close to those shown, and provide an equally good description. 

The Jfty results are given in figures [22] through [27] with pp = 1.2, and are again very satis- 
factory. We have included low-energy data, although they are not very precise and have some 
contamination from nucleon breakup. Because of the mass of the charm quark the data on 
ul/ut ratio do not provide a strong constraint, and even at Q 2 = 50 GeV 2 the ratio is still far 
from its asymptotic value. The S2 results are again to be slightly preferred overall. In the case 
of the Jf& the results obtained using the wave functions derived from the specific potentials 
do not provide a good description of the data, the cross section being about a factor of three 
too large. There are also significant differences in <Jp/a T , the ratio rising almost linearly with 
increasing Q 2 , and by Q 2 = 50 GeV 2 is a factor of two larger than the results shown in figures 
[2~4l and |27|. However the present data cannot sensibly distinguish between this result and the 
one shown. That there is a significant difference in the predictions for the J/\fr is not surprising 
given the very different wave function used in the fit compared to those obtained from specific 
potentials: recall figure ||. 

There is very little difference in the energy dependence predicted by models SI and S2, so we 
shall show only the latter. The energy dependence at fixed Q 2 is shown in figures [28] and 



for the p, in figure |30] for the and in figure [31] for the J/ip, in each case for the S2 model. 
The break in the curves is because of the different value of e used at the lower energies. The 
model succeeds well in reproducing the trends of the data, and is particularly successsful for 
the and J/ip, even at Q 2 = in the latter case. The increasing energy dependence with 
increasing Q 2 is well represented by the model, reflecting the increasing importance of the hard 
pomeron as Q 2 is increased. This is not artifically imposed on the model. It occurs naturally 
via the loop integrals involving the meson wave function and the gluon propagators. The 
model also automatically takes account of the increasing importance of the hard pomeron with 
increasing quark mass. For example, for the p at Q 2 = the soft pomeron contributes 95% 
of the amplitude at y/s = 15 GeV and 85% of the amplitude at y/s = 75 GeV. At Q 2 = 20 
GeV 2 these have become 25% of the amplitude at ^fs = 15 GeV and 10% of the amplitude 
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at y/s = 75 GeV. These proportions are very similar for the <ft with a slight increase of the 
hard component: for example at Q 2 = the soft pomeron contributes 90% of the amplitude 
at y/jj = 15 GeV and 70% of the ampitude at y/s = 90 GeV. For the Jfl! at Q 2 = the 
soft pomeron contributes 80% of the amplitude at \fs = 15 GeV and 35% at ^/s = 250 GeV. 
These results for the soft pomeron proportions at Q 2 = are comparable with those obtained 
in other phenomenological approaches 0, HTj. The same general features, a slow onset of the 



perturbative region where the hard pomeron dominates, are compatible with those of |6g, [69 
In |68|] a theoretical analysis is made of ^*p — > p°p based on the BFKL formalism, and it is 
concluded that the perturbative term does not dominate until energies and virtualities beyond 
those currently accessed at HERA. In [ESI an analysis of J /if? photoproduction in a dipole model 
clearly illustrates the mixing of perturbative and nonperturbative effects in this process. 



Finally, examples of the predicted t-dependence are shown in figures 32 to 35. The data for 



the p and <fi are unnormalised, so the theoretical curves have been renormalized to the data. 
The relative normalisation of SI and S2 is that given by the models. For the J/& the data 
and the predictions are normalised. Clearly the description is again very satisfactory, with the 
possible exception of the p where the models predict a somewhat faster decrease with t than 
is observed. 



6 Conclusions 



We have presented a model of vector meson electroproduction, and also for photoproduction in 
the case of the J/Hf. The model gives a good overall description of the data for p, (ft and J/ip 
electroproduction with only five adjustable parameters. The dynamical mechanism, with two 
adjustable parameters, is common to each and the only freedom in going from one vector meson 
to another is the parameter pp which is related to the "size" of the meson in its rest frame. 
The values of pp required for each of the p, and J/ip are 0.5 GeV 2 , 0.6 GeV 2 and 1.2 GeV 2 
respectively, in accord with what one would expect. Not surprisingly, electroproduction of the p 
is primarily non-perturbative at small Q 2 , and an important non-perturbative component is still 
present at the highest energy and largest Q 2 for which data exist. In contrast the perturbative 
contribution dominates in J/ip production at high energy, although at small Q 2 interference 
with the non-perturbative contribution remains important. For W > 50 GeV, Q 2 > 20 GeV 2 
J/ty electroproduction can be considered to be exclusively perturbative. We finally note that 
this type of model could be used to give an excellent description of any of the p, (f> and J ftp if 
they were to be considered in isolation. 
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Figure 10: Q 2 dependence of the p-meson Figure 11: Q 2 dependence of the p-meson 
cross-section at W = 75 GeV in model SI. The cross-section at W = 15 GeV in model SI. The 
data are from: HI pi; and ZEUS [M PJ O. data are from: NMC Kg . 
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Figure 13: Q 2 dependence of the p-meson Figure 14: Q 2 dependence of the p-meson 
cross-section at W = 75 GeV in model S2. The cross-section at W = 15 GeV in model S2. The 
data are from: HI pi; and ZEUS [M PJ O. data are from: NMC Kg . 
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Figure 15: Q 2 dependence of the p-meson 
longtitudinal to transverse cross-section ratio 
at W = 75 GeV in model S2. The data are 
from: CHIO gj; NMC gSJ; E665 |7[; HI 
HI; and ZEUS @ H ||. 
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Figure 18: Q 2 dependence of the -meson 
longtitudinal to transverse cross-section ratio 
at W = 90 GeV in model SI. The data are 
from: HI 0; and ZEUS gl . 
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Figure 19: Q 2 dependence of the -meson Figure 20: Q 2 dependence of the -meson 
cross-section at W = 90 GeV in model S2. cross-section at W = 14 GeV in model S2. The 
The data are from: HI 
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Figure 21: Q 2 dependence of the -meson 
longtitudinal to transverse cross-section ratio 
at W = 90 GeV in model S2. The data are 
from: HI 0; and ZEUS gl . 
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Figure 25: Q 2 dependence of the JfQ -meson Figure 26: Q 2 dependence of the Jf® -meson 

cross-section at W = 90 GeV in model S2. The cross-section at W = 14 GeV in model S2. The 
data are from HI |54j and ZEUS 



data are from The data are from EMC 103 



E401 @; E516 H; NA14 ||; and E687 |47 . 
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Figure 27: Q 2 dependence of the Jffy -meson 
longtitudinal to transverse cross-section ratio 
at W = 90 GeV in model S2. The data are 
from: HI |>| 



331; and ZEUS 58 62 . 
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Figure 28: W dependence of the p -meson Figure 29: W dependence of the p -meson 
cross-section for various Q 2 in model S2. The cross-section for various Q 2 in model S2. The 



data are from: NMC (48|; HI @]; and ZEUS data are from HI |64]. 
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Figure 30: W dependence of the -meson Figure 31: W dependence of the J/M/ meson 
cross-section for various Q 2 in model S2. The cross-section for various Q 2 in model S2. The 



data are from: NMC g§; H1 K and ZEUS data are from: EMC @5 E401 05 E516 II 
5§ 0§. NA14 g§; E687 |7|; HI gj§; and ZEUS 
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Figure 32: p meson t dependence, (Q 2 ) = 4.8 
GeV 2 , 30 < W < 140 GeV. The dashed line 
is for the SI model, the dotted line for the S2 
model. The overall normalisation is arbitrary, 
but the relative normalisation of the SI and 
S2 predictions are as given by the model. The 



data are from HI 64 
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Figure 33: p meson t dependence, (Q 2 ) = 10.9 
GeV 2 , 30 < W < 140 GeV. The dashed line is 
for the SI model, the dotted line for the S2 
model. The overall normalisation is arbitrary, 
but the relative normalisation of the SI and 
S2 predictions are as given by the model. The 



data are from HI 64 
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Figure 34: meson t dependence, (Q 2 ) = 4.5 
GeV 2 , 40 < W < 130 GeV. The dashed line 
is for the SI model, the dotted line for the S2 
model. The overall normalisation is arbitrary, 
but the relative normalisation of the SI and 
S2 predictions are as given by the model. The 



data are from HI |67| . 
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Figure 35: J/ip meson t dependence, Q 2 = 0, 
91 < W < 110 GeV. The dashed line is for the 
SI model, the dotted line for the S2 model. 
The data are from: HI 



; and ZEUS |5| 
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